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ABSTRACT 

NASA Goddard Space Flight Center (GSFC) is 
developing a new space-bome GPS receiver that can 
operate effectively in the full range of Earth orbiting 
missions from Low Earth Orbit (LEO) to geostationary 
and beyond [1], Navigator is designed to be a fully space 
flight qualified GPS receiver optimized for fast signal 
acquisition and weak signal tracking. The fast acquisition 
capabilities provide exceptional time to first fix 
performance (TTFF) with no a priori receiver state or 
GPS almanac information, even in the presence of high 
Doppler shifts present in LEO (or near perigee in highly 
eccentric orbits). The fast acquisition capability also 
makes it feasible to implement extended correlation 
intervals and therefore significantly reduce Navigator’s 
acquisition threshold. This greatly improves GPS 
observability when the receiver is above the GPS 
constellation (and satellites must be tracked from the 
opposite side of the Earth) by providing at least 10 dB of 
increased acquisition sensitivity. Fast acquisition and 
weak signal tracking algorithms have been implemented 
and validated on a hardware development board. A fully 
functional version of the receiver, employing most of the 
flight parts, with integrated navigation software is 
expected by mid 2005. An ultimate goal of this project is 
to license the Navigator design to an industry partner who 
will then market the receiver as a commercial product. 

INTRODUCTION 

GPS has found wide application for precision spacecraft 
navigation and formation flying applications in low Earth 
orbits (LEO), but recent advances in GPS receiver designs 
and signal processing capabilities now make it feasible to 
consider using GPS to provide autonomous, onboard 
navigation capabilities for geostationary (GEO) or other 
high altitude space missions. The Geostationary 
Operational Environmental Satellite (GOES) Program 
provides one example of a program that would derive 
significant benefits from a GPS receiver capable of 
autonomously determining satellite position to 100 meters 
or less in a GEO orbit; however, options are extremely 
limited for procuring a receiver that could provide both 
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the performance and reliability required to meet the 10+ 
year lifetime of the GOES satellites. 

In recent years, many organizations have contributed to 
advancing capabilities for applying GPS to high altitude 
orbits through analytical studies, flight experiments, and 
GPS receiver development efforts. TRW (now Northrop 
Grumman) developed a highly successful GPS based 
translator system that has been used to perform 
operational orbit determination for a US military GEO 
satellite [2], Alcatel and CNES jointly developed a GPS 
receiver for the geostationary STENTOR satellite that 
employed a tightly coupled navigation filter/tracking loop 
design to enable weak signal tracking [3]. Many other 
commercial manufacturers of space capable GPS 
receivers including General Dynamics and Surrey 
Satellite Technology have recently considered plans to 
market space receivers for high altitude applications. The 
Mission Engineering and Systems Analysis Division at 
Goddard Space Flight Center (GSFC) previously 
developed a GPS receiver called PiVoT (Figure 1) based 
on the commercially available GPS development system 
marketed by Zarlink (formerly Mitel and before that GEC 
Plessey). The PiVoT receiver was originally developed 
for LEO applications; however, it incorporates a number 
of software modifications that enable its operation above 
the GPS constellation, and it has been used extensively at 
GSFC for hardware in-the-loop testing in high altitude 
orbits. PiVoT uses a high quality, temperature- 
compensated crystal oscillator that is integrated with the 
GPS Enhanced Orbit Navigation System (GEONS) 
software. 



Figure 1: The PiVoT GPS Receiver 


This paper describes a project at GSFC to develop a new, 
fully space qualified GPS receiver targeted for high 
altitude applications. The Navigator receiver project 
builds on GSFC’s extensive experience in high altitude 
GPS applications, and has been developed around the 
following design specifications: 

• Build a fully space flight qualified GPS receiver 
using radiation-hardened (rad-hard) parts 

® Implement algorithms in FPGAs to allow for easy 
modification, upgrading, and customization 

• Acquire and track signals at Carrier to Noise Density 
Ratio, C/N 0 325 dB-Hz 


• Acquire GPS signals (with no a priori knowledge) 
within one second for strong signals (>40 dB-Hz) and 
one minute for weak signals (<40 dB-Hz) 

Two key performance improvements over the PiVoT 
receiver are: 1) signals can be reliably acquired and 
tracked down to a C/N 0 of approximately 25 dB-Hz, a 
10 dB improvement in sensitivity; and 2) Fast Fourier 
Transform (FFT) based acquisition algorithms provide 
extremely short acquisition times even for the weakest 
signals. Figure 2 is a picture of the Navigator 
development board. 




Figure 2: Navigator GPS Development Board 


The remainder of the paper is organized as follows. First 
some theoretical background is provided for the 
acquisition and tracking algorithms implemented in 
Navigator. Next the various functions and performance 
of the algorithms are described in detail, and the process 
of porting algorithms from software to a hardware 
development board is discussed. The actual flight 
hardware and receiver software designs are discussed, and 
finally, the current status of the project is summarized. 

MOTIVATION AND THEORETICAL 
BACKGROUND 

In order to make use of a GPS signal coming from a 
particular satellite, a GPS receiver must first acquire and 
then track it. Acquisition is by far the most 
computationally demanding task, requiring a search 
across a three dimensional space of unknown time delay, 
Doppler shift, and satellite number. It is also the limiting 
factor for sensitivity for the common GPS receiver 
architecture. Tracking of signals down to 25dB-Hz is 
possible with standard GPS methods [10], while 
acquisition is limited to about 35dB-Hz. For Navigator to 
achieve its design goals, a new acquisition scheme was 
needed. 

The acquisition algorithms implemented in Navigator are 
in large part directly adapted from the paper by Psiaki on 
acquisition of weak GPS signals [4]. In his paper, Psiaki 


2 




presents a detection statistic and hypothesis testing 
scheme that enables the acquisition of very weak GPS 
signals with no a priori knowledge. In particular, the 
algorithms assume no knowledge of the GPS data 
message bits. He describes the use of the FFT to 
efficiently calculate this detection statistic, as well as two 
important computational devices used to reduce the 
number of FFTs that need to be executed. He also 
discusses the problem of acquiring a weak signal in the 
presence of a strong signal(s), and demonstrates one 
possible way to deal with it. 

Psiaki demonstrates the effectiveness of his algorithms 
through MATLAB implementations. However, he 
concludes that it is appropriate for use in a software 
receiver in off-line processing. In one test case, where a 
21dB-Hz signal is acquired, Psiaki reports a processing 
time of 14hrs and a 240MB storage requirement running 
the algorithm in MATLAB on a 733MHz Pentium II. In 
this paper, we describe a considerable effort undertaken to 
tailor these algorithms for implementation in rad-hard 
FPGAs. This allows the algorithms to run in real-time 
and to achieve the desired sensitivity with acquisition 
times of only seconds. 

The remainder of this section reviews the theoretical 
motivation for the methods used by Navigator to achieve 
its improved sensitivity. See reference [4] as well. 

Signal Model and Test 

The received, downconverted, and sampled LI C/A GPS 
signal can, assuming only one signal in view for now, be 
modeled as: 

y k =As k +n k 

s k = c w (kT s -r)e JwkT ‘ = c k (r,o>)e J ^ (1) 

A=\A\ e je = a^T s e j \n k ~ C/V(0,cr 2 ) 

The received signal y k , consists of the sum of the GPS 
signal, A s k (with complex amplitude A related to C/N 0 as 
shown with T s the sampling period in sec.), and complex 
white Gaussian noise, n k (real and imaginary part of n k are 
N(0,a 2 /2) and independent). The noise can be modeled as 
white because we assume that the signal is bandpass 
filtered and then sampled at the Nyquist rate for the 
resulting bandwidth. s k depends on two parameters 
(which are assumed constant over k=0..N-l): co, the 
Doppler shift in rad/s and t, the time delay in sec. c ro (.) is 
the Doppler compressed/expanded version of the 
1 ,023Mchip/s, 1ms (1023 chip) periodic PRN spreading 
code c(.), which takes values in {-1,+1}. The relation 
between them is given by c m (x)=c(x*(l+ro/co L i))- We will 
write c k (x,co) = c 0 ,(kT s -x) and normally shorten it to just 
c k . 

s k is additionally modulated by the GPS data message, 
d(.), a nominally 50bits/sec, {-1,+1}, data sequence 
whose transitions are synchronized to the beginning of 


each 20 th period of c(.). The d(.) complicates the 
acquisition problem, and we will first address the problem 
pretending it is not present. 

As mentioned, we first assume there is only one known 
satellite signal potentially present. Then, GPS signal 
acquisition consists of determining whether or not this 
signal is present and if so, it’s Doppler shift and delay, 
(co,t). x only needs to be determined modulo one code 
period since c(.) is periodic. Given a finite set of samples, 
y=(yo,yi,--.yN-i)’ (s, n defined analogously), we want to 
test between the following hypotheses, and if Hi is 
accepted, estimate (co,x): 

ff 0 :y = n O MM (2) 

ff,:y = Js + n o |i|>0 

Unfortunately the signal, As, depends on a number of 
unknown parameters (A,co,x) which lie in some known 
space. Denote the true values of these parameters: 
(A*,<a*,x,). Because of these unknowns, no clear optimal 
test exists under any common criteria. One approach is to 
use a Generalized Likelihood Ratio Test (GLRT) [8]. This 
procedure will usually produce the optimal test if it exists 
and otherwise generally gives effective sub-optimal 
testing procedures. The GLRT prescribes the following 
test: 


GLR = 


max (^,r)/i(y;4®,r) 


/o(y) 


>r 


( 3 ) 


(•) denotes the joint probability density of the data 

under hypothesis Hj. That is, compute the ratio of the 
corresponding joint probability densities of the data under 
each hypothesis, with the observed data plugged in (the 
likelihoods). Then we choose H! if this ratio exceeds 
some threshold, which will be chosen to fix or bound the 
false alarm probability, P FA , using the Neyman-Pearson 
criterion. Unknown parameters are replaced by their 
Maximum Likelihood (ML) estimates under Hi. The ML 
estimate of A is easily determined, as a function of to and 
x, to be: 

A ml(®G) = iZrto yk c k^^y JakTs (4) 


The maximization of (3) over (co,x) is usually done by a 
grid search over the range of allowable values. It can be 
shown by plugging (4) back into (3), that the GLRT 
reduces to the familiar test on the correlation statistic 
(scaling term included to normalize the noise): 


I z (.,>(y) I >r 

i=0 


(5) 


That is, look for the maximum magnitude-squared, over 
(co,x), of the correlation between the input and a unit 
amplitude, zero carrier phase replica of the signal and 



compare it against a threshold. This is equivalent to 
comparing the full ML estimate of the signal amplitude 
against a threshold, which is intuitively satisfying. It is 
also true that the ML estimates of (co,x) are provided as 
the arguments of the maximization. 

Parallel vs. Serial Search 

The test described above requires the calculation of 
correlation at each (co,x) (on a grid of test points) using 
the same set of input samples. This will be referred to as 
the parallel search. In contrast, most traditional GPS 
receivers employ a serial search. In the serial case, a local 
signal generator steps through the search grid and at each 
(co,t) grid point computes a correlation with the input 
samples as the data streams in. When a correlation 
exceeds the threshold, a detection is declared at that 
particular (o>,x). In the serial search, new data are used at 
each grid point. This amounts to performing an 
independent binary test at each grid point with (co,x) 
known (thus no need to estimate them). In fact, if phase 
is modeled as random and uniformly distributed in [0, 
2n), then this is the optimal (Uniformly Most Powerful 
for |A|>0) test under the NP and Bayes criteria. In this 
case there is no maximization and the statistics of Z(to,x) 
completely specify the performance of the test. The false 
alarm probability is set for the individual tests and thus is 
actually the false alarm rate. This rate can be traded off 
against the detection probability. In the serial search, 
false alarms slow the acquisition process, but are not 
devastating because the tracking eventually fails when 
initialized with a false signal. After the failure, acquisition 
can continue where it left off. 

As mentioned before, in the parallel search, a single block 
of data is used to maximize |Z(co,x)| over the test grid. 
Although the same computations are done as in the serial 
case (only on a fixed block of data rather than new 
blocks), the problem is statistically very different. In this 
case, the test statistic is max|Z(co,x)| whose exact 
distribution is very difficult to obtain since there is 
dependence between the correlations across the search 
grid. The statistics of Z((o,x), which fully characterize the 
serial test, can be used in the parallel case to obtain 
bounds on the false alarm probability as well as the 
approximate performance. 

One (overly-conservative) method for controlling P FA , in 
the parallel search, is via the union bound. Here, the 
desired overall false alarm probability is divided by the 
number of grid points and then used to set the threshold 
assuming (<b,t) known. A less conservative method is to 
assume: 

Pr(max (0) r) | Z(®,r) |> y) » Pr(max (ffl , | Z(cd , r') |> y) (6) 

where the (co’,x’) are the mutually independent points on 
the grid. There are approximately K^p = 
1023*DopplerRange*NT s of these independent points [4]. 
With independence, we can solve for the exact threshold 


for a given P FA . Using P FA /K indep in the (ca,x) known test 
will give a very good approximation to the desired 
threshold. 

False alarms are less tolerable in the parallel case. There 
is a range of threshold values where there is a very high 
probability that at the true signal point, Z(©»,x») will 
exceed the threshold, and so the signal should eventually 
be detected by the serial search. However, there is also a 
relatively large false alarm rate and it turns out there is a 
large probability that Z(©„x„) is not the largest correlation 
over the grid. One way to look at it is that the correct 
decision, H b is made but the estimation of (cd,x) is bad. 
We refer to this as a “Type III” error, as opposed to false 
alarms under H 0 or “Type I” errors, and missed detections 
under H, or “Type II” errors. 

Another possible approach to this problem is as an M-ary 
test with H (ovt ): “signal present at (©,x)”. Similar results 
can be obtained, however, we do not formulate the 
problem in this manner. 

If the parallel case has these new problems, why not then 
just run a serial search? The reason is that the serial 
search, where new data is used in each correlation, takes 
too long for weak signals. For example, using N 
corresponding to 1ms will allow for a best case 
acquisition of signals around 35dB-Hz, (with lots of false 
alarms). In this case TTFF could be on the order of 
30min. If we want a 10 times improvement in sensitivity 
we could increase the N to 10ms (as will be shown) but in 
addition we also need to increase the fineness of the grid 
in the frequency dimension by another factor of 10 (also 
to be shown). This implies a new TTFF on the order of 
3000 min! Somehow, the search needs to be parallelized. 
Either way, a large number of correlations need to be 
computed fast! 

Resolution of the Search Grid 

Even when the Type III error does not occur, and the 
maximization over the correlation grid identifies the true 
signal, it is unlikely that (©„x») falls exactly on a search 
grid point. Thus there are always losses caused by using a 
discrete search grid. The following results from 
expanding Eq. (5): 

= ^A^c t (T)c k (T.)e-^ kT - +n : - 

k = 0 / 

R(Ar,&ay) 

= x/2 (7) 

(£)*■£ I R ( At ’ A ®) I" ”, ~ CN ( 0 > 2 ) 

where (A©, Ax)=(©-©, : x-x,) is the Doppler and delay 
estimation error, and 8 is the (unimportant) phase of 
R(A©,Ax). This R(Aco,Ax) (which approximately depends 
only on Ax) is sometimes called the ambiguity function, 
and for us it specifies the necessary fineness of the grid 
upon which the maximization of Eq. (5) is done. 
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Uncertainty or error in the estimates of to and x result in a 
reduced mean of the correlation statistic and can be 
viewed as a decrease of the effective input C/N 0 . The 
cross sections of R(.,.) are given by the following 
equations which are plotted in Figure 3. (R~0 if 

|Ax|>T chip ). 


|R(Ar,0)| 2 »(l-^) 2 


| R(0,Aa>) | 2 = 


sm(Aa>T s N / 2) 
Nsin(Aa>T s 12) 


( 8 ) 



Figure 3: Ambiguity function cross-sections. Top 
with A 03=0, bottom with Ax=0. 


The rale of thumb followed in the development of 
Navigator is to keep |Aco|<l/(4NT s ), restricting loss along 
frequency axis to 0.2dB and to keep |Ax|<l/4 chip, which 
restricts loss to 2.5dB along the code-delay axis. 

Performance of Coherent Integration 

For (cq,x)=(co„x,) and assuming the GPS data message bit 
is constant (either constant +1 or -1) over the samples of 
interest, the statistics of Z(o3,x) are readily determined to 
be: 

XwVftyNT, 

(y) = ^e jS + n z ~ CN{JIe jS ,2) 

h 0 : l z (a) ,r)(y) l 2 ~z 2 2 ( 0 ) 

Hr- \z ( a,,r)(y)\ 2 ~ zlW 

Under FI 0 the test statistic is a chi-squared random 
variable with 2 degrees of freedom and under FF it is a 2 nd 
degree non-central chi-squared with non-centrality 
parameter X, which is a product of twice the effective 
C/N 0 and the integration time, NT S . 

As mentioned, these statistics fully characterize the serial 
search and can be used to bound the false alarm 
probability for the parallel search. They also provide the 
probability that the true signal will cross the threshold 
under H b but do not account for the Type III event that 


there is a larger correlation at another (co,x) coordinate. 
When the overall false alarm probability is kept low, 
however, the Type III event is rare and the curves below 
do give the approximate performance. 

Performance curves for the parallel search with P F a 
divided down by 1023 (appropriate for search across all 
code delays and a 1kHz Doppler range) are shown in 
Figure 4. 



Figure 4: P D vs. coherent integration time (co,x) 
known and P FA =0.05/1023 

Note that in this case, a factor of 10 increase in integration 
time results in a lOdB increase in sensitivity. 

Beyond the 20ms GPS Data bit 

Unfortunately, the 20ms GPS data bit limits how long the 
coherent integration can ran (other effects such as signal 
and clock dynamics limit this time as well but on a longer 
scale than 20ms). Correlating over a data bit transition 
will result in unpredictable cancellation of correlation 
power. Even using 20ms requires knowledge of the bit 
start times. The bit start time could be added to the 
problem as an unknown parameter and one could further 
maximize Eq. (3) over the additional dimension. This 
approach is investigated by Psiaki in his “full-bits” 
method [4]. Regardless, a single 20ms coherent 
correlation will not reliably detect the desired 25dB-Hz 
signal, as can be seen in Figure 4, so we need a way to 
extend the correlations beyond the data bits. 

One way around the problem has been suggested, also in 
[4] as the “half-bits” method, is to collect consecutive 
10ms correlations. If this collection is viewed as two 
alternating sets, then one of the two sets are guaranteed to 
be free of bit transitions. 

A collection of M coherent correlations that individually 
avoid data bit transitions, {Z (o)jt)jm , m=l,2..M}, can be 
combined “non-coherently” by taking the magnitude 
squared and then summing. This removes the problematic 
data bit. 
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M - 1 

Q{a>, t) / . I <D,r),m I 

( 10 ) 

Q* = max ( o,, r) £W> 

Q* is then used as the test statistic to be compared against 
a threshold. This is the “P iong ” detection statistic Psiaki 
suggests in [4] and justifies as the Locally Most Powerful 
(LMP) test under his assumptions. This type of combined 
coherent/non-coherent integration was previously 
suggested for GPS signal acquisition by [10] and has been 
commonly used since. Throughout the rest of the paper, 
this method of non-coherently combining M coherent 
integrations of duration L ms (i.e. N=(L*0.001)/T S 
samples) each will be referred to as the “L/M integration.” 

Another important reason to keep the coherent integration 
period short is because of the structure of the ambiguity 
function along the frequency axis, as discussed above. 
Correlations can be combined non-coherently for as long 
as desired without the need for increasing frequency 
resolution. For 10ms coherent integration, a 25Hz grid 
spacing will result in a worst case loss of only 0.2dB. 

Performance of Non-coherent Integration 

Again, for (co,t) =((o„-r») and assuming data-bit transitions 
have been avoided in each coherent integration interval, 
the statistics of Q are readily determined to be: 

Ho- Q(o,t) (y) ~ % 2 M (0) ( J 

Hr- e M (y)~zL(m 

Performance curves for the parallel search with P FA 
divided down by 1023 are shown in Figure 5. 



total integration time (sec) 

Figure 5: Pd vs. total integration time for 10/M 
integrations with (co,t) known and P FA =0.05/1023 

As may be expected, the gains from non-coherent 
integration come more slowly than from coherent 
integration. Sensitivity increases roughly with the square 


root of integration time. A 10 times longer integration 
results in only about a 5dB increase in sensitivity. 

However, the above curves are very promising. They 
suggest a 25dB-Hz signal can be very reliably acquired 
(>0.9 probability) by summing the squared magnitudes of 
only 10- 10ms coherent correlations. This is only one 
tenth of a second of data. Furthermore, 20dB-Hz and 
lower are reachable with not unreasonably long data 
records. 

As a check on the reliability of these results Figure 6 
shows the results of a small simulation study which 
examines the Receiver Operating Characteristic (ROC, P D 
vs. P FA ) for a 25dB-Hz GPS signal with zero Doppler 
assumed. 



1023*P„ 

Figure 6: ROC for 10/10 integration, 25dB-Iiz signal 
co = 0 known, x unknown, P FA divided down by 1023 


The observed curve shows good agreement with the 
theoretical result. The Type III error deviations are also 
clearly seen to be a real phenomenon but as mentioned, 
they are rare when the P FA is low and thus the threshold is 
set high. 

Interfering Signals 

The previous results are obtained under the (almost 
always false) assumption that only one signal is in view. 
Tire reason the assumption was made is that ideally the 
different PRN codes have zero cross-correlation so that 
the test-statistics would give the same result with or 
without the interfering signals. This, however, implies we 
knew what form our test/estimator would take from the 
start, and thus hurts our claim of a ground up approach. 
Reference [11] takes a more realistic approach to the 
problem, however the resulting solution is much more 
computationally demanding. 

Regardless, the exercise has still been valuable because 
this is a very effective acquisition strategy and the 
performance results hold up well when there are also 
signals of equal or lesser strength present. Unfortunately, 





the performance results do not hold up when many strong 
signals are present as well. This is due to the non-ideal 
cross-correlation properties of the C/A codes. The cross- 
correlation rejection for the C/A code is known to be 
about 23dB [9], 

The primary effect of interfering signals is that the cross- 
correlation peaks cause the number of false alarms to 
increase for a given threshold level. The simplest solution 
is to increase the thresholds to control the P FA and 
correspondingly increase the integration intervals to 
recover P D . This will work when the power difference is 
moderate, but will fail for huge disparities (>15-20dB), at 
which point the acquisition of the weak signal becomes 
hopeless without somehow dealing with the interfering 
signal. 

One possible solution is through careful selection among 
the available antennas on the spacecraft or through the use 
of antenna phasing to help reject the unwanted strong 
“jamming” signal. Alternatively, one can try to acquire 
and track the strong signals first and then cancel them out 
of the input before attempting the weak signals. This is 
suggested by [9], demonstrated with success by Psiaki [4] 
and is what is effectively done in the ML GPS Receiver 
of [11]. During tracking, precise estimates of (co,t,9) are 
maintained and all that is additionally needed is the signal 
amplitude, |A|, whose ML estimate was given in Eq. (4). 
It is possible to use this in online processing and it may be 
incorporated in the Navigator design in the future. 

NAVIGATOR RECEIVER ALGORITHMS 

The Navigator receiver can quickly and reliably acquire 
and track signals down to 25dB-Hz and lower by 
employing special hardware that computes the Q -statistic 
described above. 

Navigator Design Principles 

The Navigator receiver design is based upon two central 
design principles. First, Navigator is intended to operate 
totally autonomously. That is, it does not require external 
data aiding or any other a priori information such as a 
current estimate of time, a recent GPS almanac or a 
converged navigation filter estimate of the receiver 
dynamics, etc. If such a priori information was available 
then Navigator may be able to make use of it, but it is not 
assumed. 


memory that can be very expensive for a space qualified 
part. Further, the performance is ultimately limited by the 
size of this memory. 

DFT-Based calculation of the correlation statistic 

The primary computational device employed by 
Navigator is the use of the DFT to calculate the 1ms 
correlations. When the FFT algorithm is used to compute 
these DFTs, large computational savings are achieved. 
Specifically, when an N-point DFT is computed using the 
FFT, the savings are of order N/log 2 (N) over direct 
computation. This method calculates all code delay 
correlations in !4 chip increments (for our 
implementation) using a single 1ms block of data. This is 
a well known technique used for GPS and more generally 
for Direct Sequence Spread Spectrum (DSSS) signal 
acquisition. See [12] and its associated references. 


Let x=(x 0 , X]...x N _i)’ be a signal and X=(X 0 , Xi...X N _i) be 
its DFT sequence (likewise for c and C). The “circular- 
correlation property” gives the following identity (where 
° means point by point multiplication of the vectors): 


N - 1 

Z n ~ ^ 1 "L- 
*=0 


= IDFT[Xo C*]„ 


( 12 ) 


To apply this to our problem, use a 1ms block of the 
baseband downconverted input signal (in parentheses 
below) as the x and the DFT of the code sequence as C, 
which can be computed offline, and compute (let x now 
represent an integer offset of the code between 0 and N- 
1 ): 

z M (y) = I ](y k e- j °* T ‘)c k+T 

(13) 

= £0 )c (k+r)moAN = IDFT[X oC'] t 

k=0 


This gives the 1ms correlation. The effect of Doppler 
shift on the code sequence is ignored for now and will be 
discussed below. To get a longer correlation, add L 
consecutive 1ms block correlations: 


LN-l 

Z M ( y)= 2 

k = 0 

L - 1 (N+\)l-l L - 1 

= Z x ( * + „ = Z IDFT ^, o c*] r 

/=0 k=Nl 1=0 


(14) 


Second, Navigator operates in real time. Data is double 
buffered up front in 1ms blocks and processed as it comes 
in. Some advocate the use of large up front storage to 
buffer the entire record of sampled data needed for 
acquisition. Processing can then occur with relaxed time 
constraints, only needing to be so fast that the (co,x) 
estimate will still be valid when processing finishes. 
Then the receiver could operate in a faster than real-time 
mode to catch up with the streaming input, or simply 
begin tracking on the most recent samples. The problem 
with this method is that it requires a large up-front 


Just repeat this calculation and add up the squared 
magnitudes to get the desired Q-statistic. The Navigator 
acquisition hardware implements exactly this process. 

With this method, the Q-statistic is calculated for the 
entire code dimension based on a single block of input 
data. Further, if many FFTs can be completed in 1ms then 
the search across the Doppler dimension can be done as 
well. To do so, just repeat this method for each frequency 
on the search grid at whatever frequency granularity 
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desired. In this way, the entire (co,t) grid (at least over a 
range of Doppler frequencies) can be searched in parallel. 

Psiaki [4] outlines two methods to reduce the number of 
DFT operations required, both of which are employed in 
Navigator. The first uses the well known frequency 
shifting property of the DFT. A reference (coarse) 
Doppler shift is removed from a 1ms block of the input 
signal and its DFT is computed. The DFT vectors of all 
multiples of 1kHz offset from this reference are then 
obtained by shifting the reference DFT vector. This 
technique reduces the number of forward DFTs needed 
but not the number of IDFTs. The other explains how to 
interpolate from a coarse grid of frequencies with spacing 
on the order of hundreds of Hz to the needed fine grid 
spacing of about a factor of 10 finer. This interpolation 
onto the fine grid is similar to the “Post-Correlation FFT” 
methods that are becoming popular for GPS acquisition. 
However, Navigator uses a 1ms coherent dump time, as 
opposed to the fractional ms dump normally used, and 
rather than using an FFT, which computes the Discrete 
Time Fourier Transform (DTFT) ordinates at the Fourier 
frequencies, 1/(LNT S ), only a few DTFT values are 
computed at desired frequencies. See reference [7], 

Doppler Shift on the Code 

The first equality in Eq. (14) holds when {c k } is exactly 
N-periodic (N is such that NT s =lms). The input should 
really be correlated with an appropriately Doppler 
compressed /expanded version of the code that would not 
be exactly N-periodic. There are a few ways around this. 
Psiaki suggests resampling the input in each Doppler bin 
to enforce exactly 1 PRN period in N samples. In this case 
the lms-periodic c(.) is the correct replica. Reference [7] 
suggests a delay correction applied in the frequency 
domain prior to the IDFT. 

Navigator’s approach is to ignore the issue initially and 
correct for the ill effects later. Navigator employs a re- 
mapping of the stored correlation grid. The effect of the 
code Doppler is that the apparent correlation peak will 
drift through the correlation grid in the delay dimension at 
a rate determined by the Doppler bin under test, given by 
Eq. (15). It has been found to be sufficient to correct this 
issue only upon accumulation into long term memory 
every 20ms (via re-mapping back to appropriate code bin 
when the drift exceeds Vi chip). 

Drift Rate= — ^ chips/sec (15) 

fu 

where f chip =1.023 Mchips/sec, f u =l. 57542GHz and f dopp 
is the Doppler frequency in Hz. Assuming a 8kHz 
maximum Doppler magnitude, appropriate for GEO, this 
drift is limited by 0.1039 chips in 20ms. This actually 
serves to average out the worst (and best) case loss caused 
by code misalignment for very long acquisitions. Using 
Vi chip search spacing, rather than 2.5dB worst case loss, 
the result is an average loss of 1.2dB for all long 
integrations. 


Strong Signal Mode 

Navigator offers an extremely fast strong signal mode. In 
this mode, the acquisition is based on a single 1ms block 
of input samples. This 1ms block is buffered and then 
FFT based correlations sweep across the entire Doppler 
space at the desired granularity. The frequency shifting 
device is used in this mode, but interpolation onto the fine 
grid is not needed. 

The maximum correlation value of this operation can be 
determined without the need to store the entire grid. 
Thus, the need for a large amount of high bandwidth 
memory is entirely eliminated in this mode. 

It must be mentioned that using the maximum over the 
correlation grid as the detection statistic will necessarily 
limit the sensitivity of this mode to around 40dB-Hz. The 
threshold must be set rather high to limit the Type III 
errors. Other receivers may be able to acquire signals at 
35dB-Hz based on the same 1ms correlation but they do 
this at the cost of very high false alarm rates. For 
example, the GEC Plessey GPS Builder-2 development 
board Designer’s Guide, suggests allowing a 17% false 
alarm rate during the acquisition search, achieving a 
minimum detectable signal of around 35dB-Hz. 

Weak Signal Mode 

Navigator also offers a weak signal mode. In this mode, 
10ms coherent integrations are combined (separated by a 
10 ms delay) non-coherently by summing the squared 
magnitudes as in the definition of Q in Eq. (10). This can 
be done indefinitely, limited only by signal and clock 
dynamics. In practice, M=10 will acquire the 25dB-Hz 
signal reliably and M=100 can get 20dB-Hz. 

This mode requires two separate memories. The short 
term memory stores the 10 individual 1ms correlations for 
each coarse Doppler and the long term memory stores the 
accumulation of the Q-statistic. 

Ideally, two Q-statistics would be collected, one for each 
set of alternating 10ms blocks. In practice, however, only 
one is computed at a time. The “off’ 10ms is needed for 
computation. The second set is checked only if the first 
set fails to produce a detection. 

Since the coherent integration period is 10ms, most often 
a 25Hz Doppler spacing is used, which restricts C/N 0 
losses to 0.2 dB. Due to the fine granularity, this mode 
uses both the frequency shifting property of the DFT and 
interpolation onto a fine Doppler grid to reduce the 
number of FFT operations needed [4]. 

The maximum Doppler range that can be searched at once 
is limited by the number of FFT operations that can be 
completed in 1ms and the available bandwidth to the off 
chip memories. This will be discussed more in the 
hardware section. To cover the full Doppler range, the 
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frequencies are searched sequentially in this maximum 
Doppler block size. 

Arbitrary L/M modes are available as well, but Navigator 
most commonly uses 1/1 (strong) and 10/M (weak). For 
example, if data-aiding is available then L>10 may be 
desirable. 

Fine Acquisition 

In weak signal mode, the (co,x) estimate provided by the 
acquisition module may not be accurate enough to 
initialize tracking, particularly in the co dimension. 
Furthermore, determination of the location of the GPS 
data bit transition is essential when tracking weak signals. 


When the acquisition terminates with a successful 
detection, the estimated (cd,t) is used to initialize a 
tracking correlator/channel. Open loop 1-PRN (~lms) 
period accumulations are collected for a specified amount 
of time, e.g. J-PRN periods, during which the following 
statistic is calculated to determine an estimate of the bit 
transition time, i.e. to get bit lock (K=floor(J/20)). 


r 


= arg max 

te{0,..19} 



204+19+/ 

j=20k+l 



(16) 


with correlations extended to 20ms for tracking weak 
signals. It is shown in reference [10] that PLL phase 
tracking to 25dB-Hz is achievable with 20ms correlations 
and a good oscillator. Carrier phase tracking is known to 
be the weak link as compared with carrier frequency and 
code tracking. The theoretical bit error probability at 
C/N 0 =25dB-Hz is ~10‘ 4 for (perfect) coherent 
demodulation of BPSK signals, so data demodulation can 
be done relatively well. Much below 25dB-Hz, using 
traditional FLL/PLL/DLL tracking, reliable data 
demodulation becomes difficult and phase tracking begins 
to fail. 

Bit error probabilities can be improved by averaging 
many cycles of the repeating GPS data message. This 
method may be employed in Navigator in the future. 

In [5], Psiaki demonstrates reliable zero a priori signal 
tracking and data demodulation of the C/A signal down to 
15dB-Hz using extended Kalman filter methods. These 
methods have been thoroughly investigated by the 
authors, but were deemed to be too computationally 
intensive for the first implementation of Navigator. 

NAVIGATOR HARDWARE DESIGN 


This assumes that the initial frequency estimate is good 
enough so residual frequency error does not wash out the 
above 20ms correlations. Following the aforementioned 
rule of thumb implies assuring |Aco|<12.5Hz before 
attempting to compute Eq. (16). 

To resolve frequency more finely, we can process the 
same type of {Z k } post-correlation sequence. When the 
replica phase is kept continuous over K consecutive L ms 
coherent integration blocks, the individual correlations 
have the form: 

Z, = bd k e jWNrj + n k (17) 

Here, d k is the GPS data bit and b is a complex constant 
(essentially) independent of k. This being so, if we square 
each Z k (to remove the data bit), then we get a pure 
complex exponential plus (non-Gaussian) noise, for 
which, the maximum of the periodogram is an often used 
estimator of frequency. We can use the FFT engine to 
compute the periodogram and maximize it to get an 
estimate of 2Ao). This has been shown to work nicely in 
the lab, using L = 1ms, for signals even below 25dB-Hz 
with J = 2048, giving a resolution of <lHz for Aco. For 
smaller J, we zero-pad the Z-squared sequence. This 
assumes the frequency rate is small enough so the residual 
Doppler is roughly constant over J*L-ms. 

Depending on the situation, the two calculations can be 
done simultaneously, or first (16) then (17), or vice-versa. 

Tracking 

Navigator employs standard FLL/PLL/DLL tracking 
methods using 1ms correlations for strong signals and 


The Navigator receiver employs the conventional GPS 
receiver design consisting of a bank of hardware 
correlators controlled by a general purpose 
microprocessor. In addition, it adds a specialized 
acquisition module that rapidly calculates the long term 
detection statistic, Q*. A block diagram of the top-level 
Navigator design is shown in Figure 9. 



Figure 9: System Block Diagram 


All of the GPS specific hardware is implemented in 
VHDL to target rad-hard FPGAs. The reason for 
targeting FPGAs rather than an ASIC is partially cost 
constraint, but more importantly to maintain flexibility. 
Navigator is expected to grow in the near future (e.g. L2C 
functionality is to be added) and will act as a stepping 
stone for future GPS receiver development. FPGA 
implementation offers the greatest flexibility for growth 
and modification of the design. This is particularly 
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important for space applications because it means that 
Navigator can be customized to best fit mission specific 
goals and requirements. 

Simulation Studies 

MATLAB and C software models of the above algorithms 
were used to evaluate performance and implement 
modifications. Later models were designed to emulate 
the hardware and were used as a guide for the hardware 
design. Specifically, sampling rate, ADC bit depth, and 
all other data path bit-depths, most importantly for the 
FFT, needed to be determined. Of course, for hardware 
implementation, these rates and depths need to be held as 
small as possible without sacrificing too much 
performance. 

The sampling rate was chosen to be 2.048Msps. This is 
just over the Nyquist rate for capturing the main lobe of 
the C/A code (using complex samples) and gives 2048 
samples in the 1ms nominal code period. This allows for 
the use of the efficient radix-2 FFT algorithm. This 
choice gives a !4 chip resolution for the acquisition that 
amounts to a 2.5dB worst case loss of the effective C/N 0 . 
This is indeed a significant loss and so 4.096Msps was 
also considered. However, doubling the sampling rate 
also doubles the memory size and bandwidth requirement 
which was deemed an unacceptable cost. Psiaki [4] uses 
8.192 Msps which is even nicer performance- wise, but 
not an option for an implementation in radiation-hardened 
FPGAs. 

Reference [6] performed a thorough simulation study to 
determine the optimal bit-depths for implementation of 
GPS block processing (FFT-based acquisition) methods in 
FPGAs, although it does not consider the weak signal 
case. The study concludes that 10-bits is the optimal FFT 
bit-depth. For ADC bit-depth, it concludes that 2-bits 
results in very little performance loss but suggests that 8- 
bits is preferable in the presence of interference. 
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For Navigator, ADC depth was chosen to be 8-bits in the 
acquisition module in order to maintain a large dynamic 
range. Based on our own additional studies (including the 
weak signal case) we also conclude that 10-bits is the 
optimum bit-depth for the FFT. Simulations of 
“acquisition-margin” (ratio of largest to next largest 
correlation) vs. bit-depth show that sensitivity begins to 
fall off with fewer than 10-bits, while very little is gained 
by having more. Figure 10 shows the output of one such 
simulation. 

Acquisition HW 

The acquisition hardware accepts as input: satellite 
number, Doppler range, Doppler search granularity, and 
coherent and non-coherent integration times. It then 
calculates the Q-statistic using the following sequence of 
steps. 

1. The Acquisition FPGA receives input samples from 
the ADC. These are streamed into a 2ms double 
buffer. 

2. Samples from the ready buffer enter the Doppler 
removal block which removes the nominal carrier 
and residual (coarse) Doppler under test. 

3. Next the FFT operates on 2048 samples and the 
output is stored in a 2048 point store and shift buffer. 

4. Shifted versions of this output are multiplied against 
the pre-computed, conjugated DFT of the desired 
code and then IFFT’d to complete the 1ms 
correlation. Each one sample shift corresponds to a 
1 kHz offset from the originally removed Doppler. 

Strong signal mode acquisition finishes by sweeping 
across the Doppler range by repeating 2-4 and then 
determining the maximum magnitude (and coordinates) 
on the 1ms correlation grid which it reports back to the 
microprocessor. Figure 11 shows a block diagram of 
strong signal mode. 



3.2* 1 1 1 1 1 1 

6 7 8 9 10 11 12 

bits 

Figure 10: Acquisition margin vs. FFT 
bit depth for 25-29dB-Hz signals, 10/10 
integration, 2000 runs averaged. 


Figure 1 1 : Strong Signal Acquisition Block Diagram 
Weak signal mode finishes the following way: 

5. Each 1ms vector is stored to the short-term memory 
by repeating steps 3-4 for each coarse frequency bin 
in the search range. Programmable rounding is 
applied to allow 8-bit representation. 
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6. At the end of 10ms, during the “off’ 10ms, for each 
coarse frequency, the 10 corresponding 1ms 
correlation vectors are read from the short term 
memory. The interpolation from the coarse to the fine 
grid is applied at this point to achieve the desired fine 
frequency resolution. 

7. The squared magnitudes of the correlations are 
computed and the code Doppler correction is applied 
before it is accumulated into the long-term memory. 
Rounding here is used to maintain 16-bit words in the 
long-term memory. 

At the end of the accumulations the maximum over 
thelong-term grid, Q*, is determined and reported to the 
microprocessor for threshold comparison. 

Figure 12 shows a block diagram of weak signal mode. 



Figure 12: Weak Signal Acquisition Block Diagram 


Two potential bottlenecks limit the performance of the 
acquisition module: the bandwidth to off-chip SRAM (for 
accumulation of the correlation grids) and speed of the 
FFT operation, which is ultimately determined by the size 
of the FPGA. Running this algorithm in real-time 
requires a huge bandwidth to the off chip SRAM. To 
provide for this bandwidth, a 64-bit bus connects the 
SRAM to the acquisition FPGA. Running at 66MHz, this 
provides a 528MB/s bandwidth. Weak signal mode 
requires 32.8MB/s bandwidth per 1kHz Doppler search 
block and thus can cover 16kHz at once if limited only by 
memory bandwidth. 

The second bottleneck is in running the FFTs fast enough. 
The hardware has the capability of performing 
23FFTs/ms. Assuming a 250Hz coarse Doppler spacing, 
this implies a maximum one pass Doppler coverage of 
5.7kHz. The current design uses 4 parallel butterfly 
adders to implement the FFT, limited by the expected 
utilization of the flight FPGA. If the usage estimates are 
proven to be too conservative, the number of butterfly 
adders could be increased, providing improved FFT 
speed. 


Tracking Hardware 

The tracking FPGA consists of a standard block of 
hardware correlators. However, to improve FPGA usage 
efficiency, rather than processing samples at the sampling 
rate, data are stored in a FIFO and processed by time 
shared hardware running at a much higher rate than the 
sampler. The specific implementation consists of 3 time 
shared correlator blocks that give 12 channels each, for a 
total of 36 channels. 

Receiver Flight Part Selection 

Navigator has been designed to provide the highest levels 
of reliability in the severe radiation environment present 
in high Earth orbits: 

1. All parts must be able to withstand a Total Dose 
radiation level of 100 krad. This exposure is with no 
box shielding; spot shielding is permissible. 

2. All parts must be tolerant to a 37 MeV-cm 2 /mg 
exposure, with no SEUs (single event upsets). 

3. All parts must be SEL (single event latch-up) 
immune up to 90 MeV-cm 2 /mg. 

Selecting parts that provide the required performance and 
survivability is a significant design challenge. The flight 
RF front end will be built around Peregrine 
Semiconductor’s PE8510x L1/L2 GPS front end ASIC. 
For the FPGAs, the Actel RTAX-2000 was selected. This 
is the largest rad-hard part offered by Actel as of the 
writing of this paper. The combinatorial and sequential 
logic meet Navigator’s radiation requirements. The 
RAMs on the Actel FPGA are actually considered soft 
(not immune from SEUs); however, Actel offers different 
EDAC (Error Detection and Correction) algorithms in 
their tool set to increase the data resiliency of SRAM. 
The flight SRAM will be a 4 BAE SRAM die packaged in 
MCM (Multi-Chip Module) by 3D-Plus. 

The baseline flight oscillator will be an ovenized crystal 
oscillator (OXO); however, a version of the receiver will 
be available using a high quality temperature controlled 
crystal oscillator (TCXO) for applications that wish to 
trade some power, mass, and cost savings for slightly 
reduced performance. 

The RH-CF5208 ColdFire was chosen as the primary 
microprocessor. It is a true embedded processor with 
very low-power consumption and almost no glue logic. 
The primary disadvantage of the ColdFire is its relatively 
modest computational power (~60 MIPS @ 66MHz) and 
lack of a floating-point unit. This has been partially 
compensated for by a custom implementation of a floating 
point co-processor in the FPGAs. 

The Navigator project has been able to take advantage of 
the ongoing development of the Subsystem Data Node 
(SDN), for NASA’s Solar Dynamics Observatory (SDO) 
which provides the ColdFire CPU and external interfaces 
integrated on a 6U-220 Compact PCI card, shown in 
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Figure 13. The Core SDN design places the RH-CF5208 
CPU and peripherals (including MIL-STD-1553, RS-422, 
cPCI, EEPROM, and an ADC) on half of a 6U-220 
Compact PCI card. The other half is reserved for 
application specific hardware, in this case the GPS 
receiver hardware. 



Figure 13: Core SDN, right half will host the GPS 
specific hardware in the Navigator flight design. 


software has been developed in a manner that is hardware 
independent by utilizing a pipe-wall construct to 
communicate with the low-level software functions. As a 
result, the new software can be run with both the 
Navigator hardware and the PiVoT receiver cards, which 
has allowed Navigator software and hardware 
development to proceed in parallel. This feature also 
opens the possibility of using the Navigator software with 
other hardware platforms in the future. 



Figure 14: The Navigator software and its interface to 
the hardware. 


NAVIGATOR SOFTWARE 

The development of the Navigator receiver software 
leverages Goddard’s extensive experience with the PiVoT 
receiver, but due to licensing concerns does not rely on 
any legacy code from PiVoT. The Navigator software 
consists of two major components that have been 
developed from the ground up. A block diagram of the 
major software components is shown in 
Figure 14. 

The first component of the software is the low-level 
functions that interface directly with Navigator’s unique 
hardware to produce raw code phase, carrier Doppler, and 
carrier phase measurements. These functions include the 
control of the acquisition engine as well as the tracking 
loops, which are also specialized for the weak signal 
environment. These tasks operate at the interrupt level 
and have direct communication with the hardware. 

The second component is the basic navigation software, 
which is mn on the Nucleus real-time operating system. 
This software forms measurements, provides the standard 
position, velocity, and time point solutions when four or 
more satellites are being tracked, and handles 
commanding and telemetry messages. Additionally, the 
Navigator software includes an attitude determination 
capability when setup with a suitable antenna 
configuration. The GPS Enhanced Onboard Navigation 
System (GEONS) is integrated with the receiver software 
to provide onboard orbit determination capabilities, and to 
provide accurate state estimation/propagation during 
periods of limited/no GPS observability. The navigation 


DEVELOPMENT ENVIRONMENT AND 
STATUS 

Development Environment 

The initial hardware development was done on a 
development board procured from General Dynamics, 
shown in Figure 2. The board consists of two Xilinx 
Virtex 2000E FPGAs controlled by a Motorola Coldfire 
MCF5307 microprocessor. The board has a footprint for 
General Dynamics’ radiation hardened version of the 
MCF5307, which is the target processor for the flight 
receiver. The board also has 1MB total of SRAM and 
2MB total of FLASH PROM. Communication to a local 
PC can be accomplished via two RS-232 serial ports, over 
the debugger connection or via a high speed one-way data 
link from spare FPGA output pins to a National 
Instruments DAQ card in the PC. For the development of 
the weak signal acquisition methods we have added a 
custom SRAM daughter card with 64-bit access from one 
of the FPGAs. 

The development RF front end is built from discrete 
components. A single stage downconversion from LI 
down to 35.42MHz is done by an image-reject mixer 
driven by a high quality bench top oscillator. At this IF, 
the signal is filtered by a Zarlink SAW filter to an 
approximately 2MHz bandwidth. The signal is phase 
split and under-sampled by dual ADCs running at 
2.048MHz, leaving the final signal centered nominally at 
604kHz. 
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Limiting Errors through Simulation and Debug 

One primary goal in development was to eliminate any 
errors in translating the C-coded algorithms to VHDL. In 
order to do so, a C/VHDL co-simulation environment was 
developed. This environment consisted of three main 
parts: (1) the C program, which was actively running as a 
process on our test PC, (2) the VHDL under development, 
which was running in simulation as a second process on 
the PC, and (3) a gasket (implemented as a software pipe) 
to allow the two processes to talk to one another. 

A common data file was fed into both the C and the 
VHDL processes. Through the gasket the results of the 
tracking and acquisition were compared every millisecond 
and any discrepancies were immediately flagged. Each 
simulation generated a log file describing its behavior any 
particular data transfers, any interrupts and, any errors. 

Each logic element in the acquisition FPGA contains an 
upstream and downstream debug RAM block. These 
RAMs are used to determine the proper execution of the 
algorithm. Any element (beginning with the Doppler 
Removal Block and ending with the IFFT) can be tested 
to ensure that it is behaving properly. 

Status 

Development of the Navigator receiver is continuing 
using the commercial ColdFire and Xilinx FPGAs 
implemented on the development board shown in Figure 

2. The receiver is being tested using GPS signals from 
multiple sources: a basic GPS signal generator capable of 
generating a single GPS channel and dynamics, a Spirent 
STR4500 GPS simulator, and live GPS signals from a 
rooftop antenna. 

Fast signal acquisition and tracking has been successfully 
demonstrated on strong signals (C/N 0 >40dB-Hz) with 
LEO dynamics. For a Doppler uncertainty range of 
±50kHz and using a frequency search granularity of 
250Hz, the receiver acquires satellites on 12 channels in 
less than 2 seconds with zero a priori knowledge. Weak 
signal acquisition has been shown to acquire 25dB-Hz 
signals reliably using the 10/10 integration. Considerably 
weaker signals have also been acquired using longer 
integrations. In weak signal mode, the Doppler range that 
can be covered in one pass is currently 5.7kHz at 25Hz 
granularity. 

For development, Navigator will offer 36 tracking 
channels -and utilize 60% of one Xilinx FPGA, while 
acquisition will use approximately 75% of a second 
FPGA. The rad-hard Actel FPGAs offers less logic and 
RAM then the Xilinx counterparts, and thus will have 
higher utilization rates. For the flight design, acquisition 
will use two Actel FPGAs each at 40% utilization, while 
tracking will use 70% of a third FPGA. As the sum of the 
utilizations of the two flight acquisition FPGAs do not 
total 100%, it is possible that the acquisition design may 
fit in a single FPGA. However, the design still plans for 


three Actels total, which allows room for growth or 
features. For example, a floating point coprocessor has 
been implemented in the FPGAs to assist the ColdFire in 
its computational tasks. 

The “breadboard” version of the Navigator hardware will 
employ most of the flight parts. In particular, Actels 
(commercial AX2000’s, not the rad-hard version) rather 
than Xilinx, the rad-hard ColdFire and SRAM, the RF 
front end and the OXO. Currently, the schematic capture 
for the breadboard design is underway. The breadboard 
version is expected to be ready in mid 2005. 

SUMMARY 

Navigator is being developed as a fully space flight 
qualified GPS receiver optimized for fast signal 
acquisition and weak signal tracking. Considerable effort 
has been put forth to implement what is generally 
considered a non-real time software GPS receiver 
algorithm into real time hardware. The fast acquisition 
capabilities provide exceptional TTFF performance with 
no a priori receiver state or GPS almanac information, 
even in the presence of high Doppler shifts associated 
with LEO (or near perigee in highly eccentric orbits). 
The fast acquisition capability also makes it feasible to 
implement extended correlation intervals and therefore 
significantly reduce Navigator’s acquisition threshold. 
Fast acquisition and weak signal tracking algorithms have 
been implemented and validated on a hardware 
development board. The ultimate goal is to license the 
Navigator design to an industry partner who will then 
market the receiver as a commercial product. A fully 
functional “breadboard” version of the receiver with 
integrated navigation software is expected by mid 2005. 
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